Rotational States of Magnetic Molecules 
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We study a magnetic molecule that exhibits spin tunneling and is free to rotate about its 
anisotropy axis. Exact low-energy eigenstates of the molecule that are superpositions of spin and 
rotational states are obtained. We show that parameter a — 2(hS) 2 /(IA) determines the ground 
state of the molecule. Here hS is the spin, I is the moment of inertia, and A is the tunnel splitting. 
The magnetic moment of the molecule is zero at a < ct c — [1 — 1/(2S) 2 ]~ 1 and non-zero at a > a c . 
At a — *• oo the spin of the molecule localizes in one of the directions along the anisotropy axis. 

PACS numbers: 75.50.Xx, 33.20.Sn, 85.65.+h 



Crystals of high-spin magnetic molecules came to the 
attention of physicists after Sessoli et al. [l| discovered 
that they behave as regular arrays of identical superpara- 
magnetic particles Q • The remarkable property of mag- 
netic molecules is that their spin can tunnel between up 
and down directions Q. This leads to a characteristic 
step- wise magnetization curve discovered by Friedman et 
al. in Mni2-Acetate [I| and later observed in hundreds 
of other molecular magnets. More recently, experiments 
were performed with magnetic molecules deposited on 
surfaces [5| and with single molecules bridged between 
metallic electrodes [f|. The interest to such studies has 
been driven in part by the prospect of using magnetic 
molecules as qubits [7(. At first glance, partial or to- 
tal decoupling of the molecule from the environment ap- 
pears desirable to achieve low decoherence. It was no- 
ticed [a, 13]) however, that such a decoupling may prohibit 
spin tunneling altogether due to conservation of the total 
angular momentum, J = S + L, with S being the spin of 
the molecule and L being the orbital angular momentum 
associated with the mechanical rotation. This situation 
can be relevant to recent experiments with single mag- 
netic molecules if the molecules maintain some degree of 
freedom with respect to rotations. 

Consider first a high-spin molecule in a crystal. The 
general form of the spin Hamiltonian of the molecule is 



(1) 



where Hn commutes with S z and H± is a perturbation 
that does not commute with S z . The existence of the 
magnetic anisotropy axis means that the \±S) eigenstates 
of S z are degenerate ground states of Hn . Operator H± 
slightly perturbs the | ± S) states, adding to them small 
contributions of other \ms) states. We shall call these 
degenerate normalized perturbed states \ip±s)- Physi- 
cally they describe the magnetic moment of the molecule 
looking in one of the two directions along the anisotropy 
axis. Full perturbation theory with account of the de- 
generacy of Hs provides quantum tunneling between the 
\yj±s) states. The ground state and the first excited state 



become 



They satisfy 
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with E- — E + = A being the tunnel splitting. Since the 
crystal field Hamiltonian Hs does not possess the full in- 
variance with respect to rotations, should not be the 
eigenstates of J. However, a closed system consisting of 
the spin and the crystal does possess such invariance. It 
has been demonstrated [l(J, El that conservation of the 
total angular momentum (spin + crystal) dictates en- 
tanglement of spin states with elastic twists. This effect 
contributes to spin decoherence but does not significantly 
affect the ground state energy. 

The situation changes for a free magnetic molecule. A 
high-spin molecule usually consists of hundreds of atoms, 
making its mechanical properties similar to the mechani- 
cal properties of a tiny solid body. Free magnetic clusters 
in beams have been studied in the past [121 ]. They ex- 
hibit a number of interesting phenomena some of which 
have been attributed to the interaction between spin and 
mechanical degrees of freedom. General analytical solu- 
tion for the rotational quantum levels of a rigid body does 
not exist. Spin degree of freedom further complicates the 
problem. However, as we demonstrate below, the exact 
eigenstates and exact energy levels can be obtained ana- 
lytically for the low-energy states of a magnetic molecule 
that is free to rotate about its anisotropy axis. This 
could be the case when a free molecule is in a magnetic 
field or the molecule is bridged between two leads. The 
eigenstates of such a molecule must be the eigenstates 
of J z = S z + L z . It is then clear that, unless mechan- 
ical rotations are involved, conservation of J z prohibits 
quantum tunneling of S. For, e.g., J z = the transitions 
can only occur between the states \ifrs) <8> \ttil = —S) 
and \ip-s) ® \m,L = S). These are the states in which 
the angular momentum due to spin is compensated by 
the angular momentum due to mechanical rotation. For 
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a superposition of these states to be the ground state 
of the system, the kinetic energy, (hS) 2 /(2I), associated 
with the rotation cannot significantly exceed the energy 
gain, A/2, due to spin tunneling. Otherwise the ground 
state will be \ip±s) ® \ m L = 0). For a solid particle, the 
moment of inertia 7 grows as the fifth power of the size of 
the particle. Consequently, rotational effects should be 
less important in large particles. For magnetic molecules, 
however, the rotational energy (hS) 2 / (27) in many cases 
will be large enough to cause localization in one of the 
\ip±s) spin states. Exact analytical solution of this prob- 
lem is given below. 

Since the low-energy spin states of the molecule are 
superpositions of \ip±s), it is convenient to describe such 
a two-state system by a pseudospin 1/2. Components of 
the corresponding Pauli operator ar are 



(Ty 



\ip-s)(ips \ - 
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The projection of Hs onto \ip±s) states is 
H a = J2 (m\H s \n)\m){n\ 

m,n—tp±s 



(5) 



Expressing \if>±s) via ^± according to Eq. |2]), it is easy 
to see from Eq. ([3]) that 

(il>±s\H 3 \1>±s) = 0, (ip-s\Hs\i>s) = -A/2 . (6) 
With the help of these relations one obtains from Eq. ([5|) 
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So far we have not considered mechanical rotations of 
the molecule. Rotation by angle <fi about the anisotropy 
axis Z, transforms the spin Hamiltonian into 
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Noticing that 

S z \i>±s) 



H s e l 



S z \±S) = ±S\ij±s) 



(8) 
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it is easy to project Hamiltonian ([8]) onto ip±s- Simple 
calculation yields the following generalization of Eq. ([7]): 
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The full Hamiltonian of the magnetic molecule rotating 
about its anisotropy axis is 
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~a x cos(2SV>) + <j y sin(25</>)] , (11) 



with L z — —i(d/d(f>). 

We are now in a position to find the eigenstates of the 
rotating molecule. By construction, the Hamiltonian (fTTj) 
is invariant with respect to rotations about the Z-axis. 
Consequently, its eigenstates must be the eigenstates of 
J z = L z -\- S z - 
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S)i + Cs\il>-s)<Z>\J + S)i) ■ 

(12) 



Here J = m,j while index I denotes states in the mechan- 
ical space, with \m)i = \mi) — exp(im^0). Solution 
of H"i> j = Ej^j gives the following expression for the 
energy levels: 
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where 
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For J ^ each state is degenerate with respect to the 
sign of J. For J = 0, 1, 2, ... coefficients C± are given by 
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c^ s = tJi 
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VS 2 + {aJ) 2 
aJ 

~ ^S 2 + {aJf 



(15) 



where =F correlates with ± in Eq. p3)l . 

At J k mi > S, Eq. (fT3|) gives the energy of the me- 
chanical rotation, (hm^) 2 / (21). At small a the ground 
state and the first excited state correspond to J = 0, 
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Here the first term is the energy of the rotation with 
m,L = ±5. For a molecule rigidly coupled to an infinite 
mass one has 7 — * oo and the energy of the rotation 
goes to zero. In this case one recovers from Eq. (fT6|) 
the energies, ±A/2, of the tunnel-split spin states in a 
macroscopic crystal. As a increases, the ground state 
switches to higher J. The value of a at which the ground 
state changes from Ej_i to Ej satisfies 



Ej-x-(otj) = Ej-(atj) . 
Solution of this equation for J = 1, 2, S gives 



otj = 



1 - 
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For a smaller then 
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ai =[l- 1/(25) 
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the ground state corresponds to J = 0, C±s = 1. At 
a = a.\ the transition to the J = 1 ground state takes 
place. At a = a 2 the ground state changes from J = 1 
to J = 2, and so on. At a greater than 
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(20) 



the ground state always corresponds to J = 5. For, e.g., 
S 1 = 10 one obtains ai = 1.0025 and ai = 3.2066. The 
dependence of the ground state energy on a for S = 10 
is shown in Fig. [TJ While this dependence is smooth, 
the derivative of the ground state energy on a shows 
steps at the critical values of a given by Eq. (fl~8]) . In 
the limit of a — > the ground state energy is —A/2. 
This is the gain in energy due to spin tunneling between 
\ip±s) states in an infinitely heavy particle. In the limit of 
a ^> 1 (light particle) J = S and according to Eq. (TT3|) 
the ground state energy approaches zero as —A/ (4a). 
This corresponds to the gradual localization of the spin 
in one of the \ip±s) states. 
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FIG. 1: Zero- field ground state energy as a function of a. 
Inset shows the derivative of the ground state energy on a. 



To compute the magnetic moment of the molecule we 
notice that L z in our formalism describes the mechanical 
rotation of the molecule as a whole, not the orbital states 
of the electrons. Consequently, the magnetic moment of 
a free magnetic molecule should be entirely due to its 
spin: 



/i = -gfiB{^j\S z \^j) = -q^bS 



aj 



V$ 2 + W 



(21) 



Here g is the spin gyromagnetic factor. The minus sign 
reflects negative gyromagnetic ratio, 7 = —gfis/% for 
the electron spin. If the ground state corresponds to J = 
0, spin-up and spin-down states contribute equally to the 
wave function and the magnetic moment is zero. When 
J in the ground-state is non-zero, spin-up and spin-down 



states contribute with different weights and the molecule 
has a non-zero magnetic moment. Which J corresponds 
to the ground state depends on the parameter a. The 
dependence of the ground state magnetic moment on a 
is shown in Fig. O 




FIG. 2: Dependence of the ground-state magnetic moment on 
parameter a. 



The above results can be easily generalized to take into 
account the effect of the external magnetic field B ap- 
plied along the Z-axis. Such a field adds a Zeeman term, 
9^bS z B, to the Hamiltonian ([TJ) . This term is invariant 
with respect to the rotation by the angle <f>. Its projec- 
tion onto ip±s simply adds g^,BSBa z to Eq. (fTTj) . The 
full projected Hamiltonian becomes 



H = 



' 21 d<b 2 



where 



[a x cos(2S(f>) 



W = -2gfi B SB . 



i(2S$] 



(22) 
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Since this Hamiltonian is invariant with respect to rota- 
tions about the Z-axis, its eigenfunctions are still given 
by Eq. JT2]) with the coefficients C±s now depending on 
B. Solving H^ j = Ej^fj one obtains 



Ej± = 




(24) 



for the energy levels. Here W can be positive or negative 
depending on the orientation of the field. Positive W 
corresponds to the magnetic field in the direction of the 
magnetic moment, which provides the lower energy. At 
coefficients C±s can be presented in the form 



C s 
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with 



in the ballpark of 10 35 g cm 2 . 



W = W + aJA/S . 



(26) 



Notice that Eq. (|25f coincides with the form of C±s for 
a frozen magnetic molecule in the magnetic field B = 
B — hJ/(-fI). The magnetic moment of the molecule 
that is free to rotate is given by 
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VA 2 + W 2 



(27) 



In the absence of the magnetic field, quantum number 
J corresponding to the ground state is determined by a. 
For a given magnetic molecule this parameter is fixed. 
On the contrary, in the presence of the field, J can be 
manipulated by changing B. Solving Eq. (flTI) with Ej± 
of Eq (|2"4"]) . one obtains the following expression for W = 
Wj at which the ground state switches from J — 1 to J: 
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Here J — 1, 2, S. The dependence of the ground state 
magnetic moment on W is shown in Fig. [3] The jumps 
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FIG. 3: Field dependence of the magnetic moment. Note the 
jumps at W = Wj. 



at critical values of the field must show as sharp maxima 
in the differential susceptibility. 

In conclusion, we have obtained exact low-energy 
quantum states of a magnetic particle (molecule) that 
exhibits spin tunneling and is free to rotate about its 
anisotropy axis. The ground state depends on the param- 
eter a — 2(hS) 2 /(IA). Various limits studied above are 
physically accessible in magnetic molecules and atomic 
clusters. At a — > oo the spin localizes in one of the two 
directions along the magnetic anisotropy axis. Magnetic 
molecule of a nanometer size has the moment of inertia 



For S — 10 this provides 
a ~ 1 at A/h ~ 10 10 s _1 . The tunnel splitting of, e.g., 
Mn 12 and Feg molecules is much smaller. Thus the spin 
tunneling in these molecules must be strongly suppressed 
if they are free to rotate. This effect may be important 
in designing qubits based upon magnetic molecules. 
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